Genetic and cultural adaptations underlie the establishment of dairy pastoralism in the Tibetan Plateau.

BACKGROUND: Domestication and introduction of dairy animals facilitated the permanent human occupation of the Tibetan Plateau. Yet the history of dairy pastoralism in the Tibetan Plateau remains poorly understood. Little is known how Tibetans adapted to milk and dairy products. RESULTS: We integrated archeological evidence and genetic analysis to show the picture that the dairy ruminants, together with dogs, were introduced from West Eurasia into the Tibetan Plateau since ~ 3600 years ago. The genetic admixture between the exotic and indigenous dogs enriched the candidate lactase persistence (LP) allele 10974A > G of West Eurasian origin in Tibetan dogs. In vitro experiments demonstrate that - 13838G > A functions as a LP allele in Tibetans. Unlike multiple LP alleles presenting selective signatures in West Eurasians and South Asians, the de novo origin of Tibetan-specific LP allele - 13838G > A with low frequency (~ 6-7%) and absence of selection corresponds - 13910C > T in pastoralists across eastern Eurasia steppe. CONCLUSIONS: Results depict a novel scenario of genetic and cultural adaptations to diet and expand current understanding of the establishment of dairy pastoralism in the Tibetan Plateau.

A new snake of the genus Dendrelaphis Boulenger, 1890 (Squamata: Colubridae) from the coastal area of southern Vietnam.

We describe a new species of Dendrelaphis from the coastal area in southern Vietnam based on morphological data and nucleotide sequences from one partial mitochondrial gene (Cytb). Dendrelaphis binhi sp. nov. is diagnosed by the following morphological characters: body scale rows 13 at neck and midbody, 9-11 before vent; vertebral scale row feebly enlarged; ventrals 154-158 in males and 161-170 in females; subcaudals 95-106 in males and 95-102 in females; one supralabial touching the eye; 18 or 19 maxillary teeth; hemipenis spinose, reaching 10th or 11th SC, with a slender papilla; cloacal plate divided; dorsum dark brown; and a white stripe along the lower flank. The new species differs from its congeners by an uncorrected p-distance in Cytb sequences of at least 8.1%.

Association filtering and generative adversarial networks for predicting lncRNA-associated disease.

BACKGROUND: Long non-coding RNA (lncRNA) closely associates with numerous biological processes, and with many diseases. Therefore, lncRNA-disease association prediction helps obtain relevant biological information and understand pathogenesis, and thus better diagnose preventable diseases. RESULTS: Herein, we offer the LDAF_GAN method for predicting lncRNA-associated disease based on association filtering and generative adversarial networks. Experimentation used two types of data: lncRNA-disease associated data without lncRNA sequence features, and fused lncRNA sequence features. LDAF_GAN uses a generator and discriminator, and differs from the original GAN by the addition of a filtering operation and negative sampling. Filtering allows the generator output to filter out unassociated diseases before being fed into the discriminator. Thus, the results generated by the model focuses only on lncRNAs associated with disease. Negative sampling takes a portion of disease terms with 0 from the association matrix as negative samples, which are assumed to be unassociated with lncRNA. A regular term is added to the loss function to avoid producing a vector with all values of 1, which can fool the discriminator. Thus, the model requires that generated positive samples are close to 1, and negative samples are close to 0. The model achieved a superior fitting effect; LDAF_GAN had superior performance in predicting fivefold cross-validations on the two datasets with AUC values of 0.9265 and 0.9278, respectively. In the case study, LDAF_GAN predicted disease association for six lncRNAs-H19, MALAT1, XIST, ZFAS1, UCA1, and ZEB1-AS1-and with the top ten predictions of 100%, 80%, 90%, 90%, 100%, and 90%, respectively, which were reported by previous studies. CONCLUSION: LDAF_GAN efficiently predicts the potential association of existing lncRNAs and the potential association of new lncRNAs with diseases. The results of fivefold cross-validation, tenfold cross-validation, and case studies suggest that the model has great predictive potential for lncRNA-disease association prediction.

Revalidation and expanded description of Mustela aistoodonnivalis (Mustelidae: Carnivora) based on a multigene phylogeny and morphology.

The lacked-teeth pygmy weasel, Mustela aistoodonnivalis Wu & Kao, 1991, was originally described as being from Taibai Mountain and Zhashui county, Shaanxi, China. Subsequently, it was considered a subspecies or synonym of Mustela nivalis. In a faunal survey of northwestern Sichuan, eight specimens of M. aistoodonnivalis were collected. A molecular phylogenetic analysis of one mitochondrial and six nuclear genes clustered the specimens as a distinct clade and not with M. nivalis. Morphologically, the lack of the second lower molar differentiated them from M. nivalis, and genetic distances were typical of discrete species. These analyses confirmed that M. aistoodonnivalis is an independent species in the genus Mustela.

Population structure and evolutionary history of the greater cane rat (Thryonomys swinderianus) from the Guinean Forests of West Africa.

Grasscutter (Thryonomys swinderianus) is a large-body old world rodent found in sub-Saharan Africa. The body size and the unique taste of the meat of this major crop pest have made it a target of intense hunting and a potential consideration as a micro-livestock. However, there is insufficient knowledge on the genetic diversity of its populations across African Guinean forests. Herein, we investigated the genetic diversity, population structures and evolutionary history of seven Nigerian wild grasscutter populations together with individuals from Cameroon, Republic of Benin, and Ghana, using five mitochondrial fragments, including D-loop and cytochrome b (CYTB). D-loop haplotype diversity ranged from 0.571 (± 0.149) in Republic of Benin to 0.921 (± 0.013) in Ghana. Within Nigeria, the haplotype diversity ranged from 0.659 (± 0.059) in Cross River to 0.837 (± 0.075) in Ondo subpopulation. The fixation index (FST), haplotype frequency distribution and analysis of molecular variance revealed varying levels of population structures across populations. No significant signature of population contraction was detected in the grasscutter populations. Evolutionary analyses of CYTB suggests that South African population might have diverged from other populations about 6.1 (2.6-10.18, 95% CI) MYA. Taken together, this study reveals the population status and evolutionary history of grasscutter populations in the region.

DNA barcoding of Chinese snakes reveals hidden diversity and conservation needs.

DNA barcoding has greatly facilitated studies of taxonomy, biodiversity, biological conservation, and ecology. Here, we establish a reliable DNA barcoding library for Chinese snakes, unveiling hidden diversity with implications for taxonomy, and provide a standardized tool for conservation management. Our comprehensive study includes 1638 cytochrome c oxidase subunit I (COI) sequences from Chinese snakes that correspond to 17 families, 65 genera, 228 named species (80.6% of named species) and 36 candidate species. A barcode gap analysis reveals gaps, where all nearest neighbour distances exceed maximum intraspecific distances, in 217 named species and all candidate species. Three species-delimitation methods (ABGD, sGMYC, and sPTP) recover 320 operational taxonomic units (OTUs), of which 192 OTUs correspond to named and candidate species. Twenty-eight other named species share OTUs, such as Azemiops feae and A. kharini, Gloydius halys, G. shedaoensis, and G. intermedius, and Bungarus multicinctus and B. candidus, representing inconsistencies most probably caused by imperfect taxonomy, recent and rapid speciation, weak taxonomic signal, introgressive hybridization, and/or inadequate phylogenetic signal. In contrast, 43 species and candidate species assign to two or more OTUs due to having large intraspecific distances. If most OTUs detected in this study reflect valid species, including the 36 candidate species, then 30% more species would exist than are currently recognized. Several OTU divergences associate with known biogeographic barriers, such as the Taiwan Strait. In addition to facilitating future studies, this reliable and relatively comprehensive reference database will play an important role in the future monitoring, conservation, and management of Chinese snakes.

Divergent sensory and immune gene evolution in sea turtles with contrasting demographic and life histories.

Sea turtles represent an ancient lineage of marine vertebrates that evolved from terrestrial ancestors over 100 Mya. The genomic basis of the unique physiological and ecological traits enabling these species to thrive in diverse marine habitats remains largely unknown. Additionally, many populations have drastically declined due to anthropogenic activities over the past two centuries, and their recovery is a high global conservation priority. We generated and analyzed high-quality reference genomes for the leatherback (Dermochelys coriacea) and green (Chelonia mydas) turtles, representing the two extant sea turtle families. These genomes are highly syntenic and homologous, but localized regions of noncollinearity were associated with higher copy numbers of immune, zinc-finger, and olfactory receptor (OR) genes in green turtles, with ORs related to waterborne odorants greatly expanded in green turtles. Our findings suggest that divergent evolution of these key gene families may underlie immunological and sensory adaptations assisting navigation, occupancy of neritic versus pelagic environments, and diet specialization. Reduced collinearity was especially prevalent in microchromosomes, with greater gene content, heterozygosity, and genetic distances between species, supporting their critical role in vertebrate evolutionary adaptation. Finally, diversity and demographic histories starkly contrasted between species, indicating that leatherback turtles have had a low yet stable effective population size, exhibit extremely low diversity compared with other reptiles, and harbor a higher genetic load compared with green turtles, reinforcing concern over their persistence under future climate scenarios. These genomes provide invaluable resources for advancing our understanding of evolution and conservation best practices in an imperiled vertebrate lineage.

Large-Scale Chromosomal Changes Lead to Genome-Level Expression Alterations, Environmental Adaptation, and Speciation in the Gayal (Bos frontalis).

Determining the functional consequences of karyotypic changes is invariably challenging because evolution tends to obscure many of its own footprints, such as accumulated mutations, recombination events, and demographic perturbations. Here, we describe the assembly of a chromosome-level reference genome of the gayal (Bos frontalis) thereby revealing the structure, at base-pair-level resolution, of a telo/acrocentric-to-telo/acrocentric Robertsonian translocation (2;28) (T/A-to-T/A rob[2;28]). The absence of any reduction in the recombination rate or genetic introgression within the fusion region of gayal served to challenge the long-standing view of a role for fusion-induced meiotic dysfunction in speciation. The disproportionate increase noted in the distant interactions across pro-chr2 and pro-chr28, and the change in open-chromatin accessibility following rob(2;28), may, however, have led to the various gene expression irregularities observed in the gayal. Indeed, we found that many muscle-related genes, located synthetically on pro-chr2 and pro-chr28, exhibited significant changes in expression. This, combined with genome-scale structural variants and expression alterations in genes involved in myofibril composition, may have driven the rapid sarcomere adaptation of gayal to its rugged mountain habitat. Our findings not only suggest that large-scale chromosomal changes can lead to alterations in genome-level expression, thereby promoting both adaptation and speciation, but also illuminate novel avenues for studying the relationship between karyotype evolution and speciation.

Evolution and diversification of Mountain voles (Rodentia: Cricetidae).

The systematics of the Cricetid genus Neodon have long been fraught with uncertainty due to sampling issues and a lack of comprehensive datasets. To gain better insights into the phylogeny and evolution of Neodon, we systematically sampled Neodon across the Hengduan and Himalayan Mountains, which cover most of its range in China. Analyses of skulls, teeth, and bacular structures revealed 15 distinct patterns corresponding to 15 species of Neodon. In addition to morphological analyses, we generated a high-quality reference genome for the mountain vole and generated whole-genome sequencing data for 47 samples. Phylogenomic analyses supported the recognition of six new species, revealing a long-term underestimation of Neodon diversity. We further identified positively selected genes potentially related to high-elevation adaptation. Together, our results illuminate how climate change caused the plateau to become the centre of Neodon origin and diversification and how mountain voles have adapted to the hypoxic high-altitude plateau environment.

The highest-elevation frog provides insights into mechanisms and evolution of defenses against high UV radiation.

Defense against ultraviolet (UV) radiation exposure is essential for survival, especially in high-elevation species. Although some specific genes involved in UV response have been reported, the full view of UV defense mechanisms remains largely unexplored. Herein, we used integrated approaches to analyze UV responses in the highest-elevation frog, Nanorana parkeri. We show less damage and more efficient antioxidant activity in skin of this frog than those of its lower-elevation relatives after UV exposure. We also reveal genes related to UV defense and a corresponding temporal expression pattern in N. parkeri. Genomic and metabolomic analysis along with large-scale transcriptomic profiling revealed a time-dependent coordinated defense mechanism in N. parkeri. We also identified several microRNAs that play important regulatory roles, especially in decreasing the expression levels of cell cycle genes. Moreover, multiple defense genes (i.e., TYR for melanogenesis) exhibit positive selection with function-enhancing substitutions. Thus, both expression shifts and gene mutations contribute to UV adaptation in N. parkeri. Our work demonstrates a genetic framework for evolution of UV defense in a natural environment.

Benchmarking ultra-high molecular weight DNA preservation methods for long-read and long-range sequencing.

BACKGROUND: Studies in vertebrate genomics require sampling from a broad range of tissue types, taxa, and localities. Recent advancements in long-read and long-range genome sequencing have made it possible to produce high-quality chromosome-level genome assemblies for almost any organism. However, adequate tissue preservation for the requisite ultra-high molecular weight DNA (uHMW DNA) remains a major challenge. Here we present a comparative study of preservation methods for field and laboratory tissue sampling, across vertebrate classes and different tissue types. RESULTS: We find that storage temperature was the strongest predictor of uHMW fragment lengths. While immediate flash-freezing remains the sample preservation gold standard, samples preserved in 95% EtOH or 20-25% DMSO-EDTA showed little fragment length degradation when stored at 4°C for 6 hours. Samples in 95% EtOH or 20-25% DMSO-EDTA kept at 4°C for 1 week after dissection still yielded adequate amounts of uHMW DNA for most applications. Tissue type was a significant predictor of total DNA yield but not fragment length. Preservation solution had a smaller but significant influence on both fragment length and DNA yield. CONCLUSION: We provide sample preservation guidelines that ensure sufficient DNA integrity and amount required for use with long-read and long-range sequencing technologies across vertebrates. Our best practices generated the uHMW DNA needed for the high-quality reference genomes for phase 1 of the Vertebrate Genomes Project, whose ultimate mission is to generate chromosome-level reference genome assemblies of all ∼70,000 extant vertebrate species.

Comparative multi-locus assessment of modern Asian newts ( Cynops, Paramesotriton, and Pachytriton: Salamandridae) in southern China suggests a shared biogeographic history.

Evolutionary biologists are always interested in deciphering the geographic context of diversification, therefore they introduced the concept of comparative phylogeography, which helps to identify common mechanisms that contribute to shared genetic structures among organisms from the same region. Here, we used multi-locus genetic data along with environmental data to investigate shared phylogeographic patterns among three Asian-endemic newt genera, Cynops, Paramesotriton and Pachytriton, which occurred in montane/submontane streams or ponds in southern China. Our 222 samples from 78 localities covered the entire range of the three genera and represented the largest dataset of this group to date. We reconstructed matrilineal genealogies from two protein-coding, mitochondrial genes, and gene network from two nuclear genes. We also estimated divergence times of major cladogenetic events and used occurrence data to evaluate niche difference and similarity between lineages. Our results revealed a common basal split in all three genera that corresponds to the separation of two geographic terrains of southern China. Those ancient divergence occurred during middle to late Miocene and likely correlate with paleoclimatic fluctuations caused by the uplift of the Qinghai-Xizang (Tibet) Plateau (QTP). Particularly, the strengthening and weakening of Asian summer monsoons during the Miocene may have profoundly impacted southern China and led to repeatedly vicariance in those newts. However, despite differences in realized niches between lineages, there is no evidence for divergence of fundamental niches. Preservation of old newt matriline lineages in mountains of southern China suggests that the region acts as both museums and cradles of speciation. Based on those results, we advocate a multi-pronged protection strategy for newts in the three genera.

Discovery of a wild, genetically pure Chinese giant salamander creates new conservation opportunities.

Effective conservation of threatened biota relies on accurate assessments and scientific guidance. As an unfortunate example, Chinese giant salamanders ( Andrias, CGS) remain critically endangered in nature. Misguided conservation efforts, e.g., commercial propagation and releasing of millions of likely non-indigenous or interspecific hybrids, have further compromised conservation initiatives. Limited information on wild populations of CGS poses a significant conservation challenge. Following 18-month long field monitoring, we now report the discovery of a wild population of CGS in a closed nature reserve in Jiangxi Province, China. Genomic assessments reveal its genetic distinctiveness and do not detect genetic admixture with other species. Based on morphological and molecular evidences, we describe this CGS as a new species Andrias jiangxiensis sp. nov. This is the only known species of CGS today with a genetically pure, reproducing, in situ population. This discovery emphasizes the important role that closed nature reserves play in protecting species, and the necessity of integrating long-term field monitoring and genetic assessments. It sets a new pathway for discovering and conserving endangered species, especially for those biotas that are similarly being extirpated by anthropogenic translocations and overexploitation.

Invasive Bullfrogs Maintain MHC Polymorphism Including Alleles Associated with Chytrid Fungal Infection.

Maintenance of genetic diversity at adaptive loci may facilitate invasions by non-native species by allowing populations to adapt to novel environments, despite the loss of diversity at neutral loci that typically occurs during founder events. To evaluate this prediction, we compared genetic diversity at major histocompatibility complex (MHC) and cytochrome b (cytb) loci from 20 populations of the American bullfrog (Rana catesbeiana) across theinvasive and native ranges in North America and quantified the presence of the pathogen Batrachochytrium dendrobatidis (Bd). Compared to native populations, invasive populations had significantly higher Bd prevalence and intensity, significantly higher pairwise MHC and cytb FST, and significantly lower cytb diversity, but maintained similar levels of MHC diversity. The two most common MHC alleles (LiCA_B and Rapi_33) were associated with a significant decreased risk of Bd infection, and we detected positive selection acting on four peptide binding residues. Phylogenetic analysis suggested invasive populations likely arose from a single founding population in the American Midwest with a possible subsequent invasion in the northwest. Overall, our study suggests that the maintenance of diversity at adaptive loci may contribute to invasion success and highlights the importance of quantifying diversity at functional loci to assess the evolutionary potential of invasive populations.

A new species of Tamiops (Rodentia, Sciuridae) from Sichuan, China.

During a small mammal survey in Wanglang National Natural Reserve, Sichuan, China in AugustOctober 2018, we collected four specimens of an Asiatic striped squirrel (genus Tamiops) that could not be assigned morphologically to any described species. Phylogenetic analyses of CYTB and nuDNA sequences strongly supported them in the genus Tamiops. The p-distances of CYTB between these individuals and other species of Tamiops ranged from 17.4 to 19.1%. Morphologically, they differed from their congeners, including T. maritimus, T. mcclellandii, T. rodolphii, and T. swinhoei by pelage, skull, and baculum characteristics. Based on the results of morphological comparisons and molecular analyses, we describe the four specimens as the new species Tamiops minshanica sp. nov.

Species Persistence with Hybridization in Toad-Headed Lizards Driven by Divergent Selection and Low Recombination.

Speciation plays a central role in evolutionary studies, and particularly how reproductive isolation (RI) evolves. The origins and persistence of RI are distinct processes that require separate evaluations. Treating them separately clarifies the drivers of speciation and then it is possible to link the processes to understand large-scale patterns of diversity. Recent genomic studies have focused predominantly on how species or RI originate. However, we know little about how species persist in face of gene flow. Here, we evaluate a contact zone of two closely related toad-headed lizards (Phrynocephalus) using a chromosome-level genome assembly and population genomics. To some extent, recent asymmetric introgression from Phrynocephalus putjatai to P. vlangalii reduces their genomic differences. However, their highly divergent regions (HDRs) have heterogeneous distributions across the genomes. Functional gene annotation indicates that many genes within HDRs are involved in reproduction and RI. Compared with allopatric populations, contact areas exhibit recent divergent selection on the HDRs and a lower population recombination rate. Taken together, this implies that divergent selection and low genetic recombination help maintain RI. This study provides insights into the genomic mechanisms that drive RI and two species persistence in the face of gene flow during the late stage of speciation.

A new species of the genus Oligodon Fitzinger, 1826 (Squamata: Colubridae) from Langbian Plateau, Vietnam.

We describe a new species of kukri snake from Da Lat City and its vicinity in Lam Dong Province, southern Vietnam on the basis of morphological characters of four specimens, including two males and two females. Oligodon tuani sp. nov. differs from other congeners by a combination of the following characters: large size in adults (TL ≤ 888 mm); 19 dorsal scale rows at midbody; ventrals 173-179 in males and 187-193 in females; subcaudals 58 or 59 in males and 44 or 45 in females; presubocular present; 8 supralabials, fourth and fifth entering orbit; cloacal plate undivided; 10 maxillary teeth, posterior three enlarged; vertebral stripe present; temporal streak present but faint and interrupted; long and deeply forked hemipenes, extending to 25th subcaudal, without spines and papillae; and 13-15+3-4 dorsal blotches.

Phenotypic and genomic adaptations to the extremely high elevation in plateau zokor (Myospalax baileyi).

The evolutionary outcomes of high elevation adaptation have been extensively described. However, whether widely distributed high elevation endemic animals adopt uniform mechanisms during adaptation to different elevational environments remains unknown, especially with respect to extreme high elevation environments. To explore this, we analysed the phenotypic and genomic data of seven populations of plateau zokor (Myospalax baileyi) along elevations ranging from 2,700 to 4,300 m. Based on whole-genome sequencing data and demographic reconstruction of the evolutionary history, we show that two populations of plateau zokor living at elevations exceeding 3,700 m diverged from other populations nearly 10,000 years ago. Further, phenotypic comparisons reveal stress-dependent adaptation, as two populations living at elevations exceeding 3,700 m have elevated ratios of heart mass to body mass relative to other populations, and the highest population (4,300 m) displays alterations in erythrocytes. Correspondingly, genomic analysis of selective sweeps indicates that positive selection might contribute to the observed phenotypic alterations in these two extremely high elevation populations, with the adaptive cardiovascular phenotypes of both populations possibly evolving under the functional constrains of their common ancestral population. Taken together, phenotypic and genomic evidence demonstrates that heterogeneous stressors impact adaptations to extreme elevations and reveals stress-dependent and genetically constrained adaptation to hypoxia, collectively providing new insights into the high elevation adaptation.

Multi-locus genetic analyses of Quasipaa from throughout its distribution.

Montane frogs of the genus Quasipaa Dubois, 1992 occur from southern China to Southeast Asia (Frost 2021). Analyses of mtDNA (Cytb) and nuDNA data (Rag1, Rag2, Rhod, Tyr) for samples from 93 localities throughout its distribution yield a phylogeny. Clades A and B occur in Southeast Asia, clade C in northern Yangtze River, China, clade D in southwestern China, and clades E and F in southeastern China. Results place Q. yei within monophyletic Quasipaa and identify two new species. Based on nuDNA data, the basal split of clade A and B indicates an Indochinese origin of Quasipaa. The west-east diversification of five species across South China (Q. spinosa, Q. exilispinosa, Q. jiulongensis, Q. shini, Q. boulengeri) corresponds to topographic terrains II and III of China. Divergence of species from southeastern China (Q. shini, Q. jiulongensis, Q. spinosa, Q. exilispinosa) and southwestern China (Q. boulengeri) dates to 15.30-16.56 Ma (million years ago). A principal component analysis (PCA) and t-test involving 19 bioclimatic variables identifies significantly different environmental conditions between the two regions. Species' distribution models (SDM) for Q. spinosa and Q. boulengeri identify the best areas to be eastern and western South China, respectively. Thus, environmental variation appears to have influenced the genetic divergence and distributions of Quasipaa in South China. Mito-nuclear discordance indicates that some individuals of Q. exilispinosa and Q. spinosa hybridized historically.